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A mini subcritical nuclear reactor was designed using Monte Carlo methods. The reactor has light water as 
moderator, natural uranium as fuel, and a 239 PuBe neutron source. In the design uranium fuel was modeled in 
an arrangement of concentric rings: 8.5, 14.5,20.5 26.5, 32.5 cm-inner radius, 3 cm-thick, and 36 cm-high. 
Different models were made from a single ring of natural uranium to five rings. For each case, the neutron 
spectra, the neutron fluence distribution, the effective multiplication factor, the amplification factor, and the 
reactor power were estimated. The ambient dose equivalent rate outside the mini reactor was also estimated. 
The maximum value for the k eff (0.78) was obtained when five rings of fuel were used; this value is close 
to 0.86 which belongs to a Nuclear Chicago subcritical reactor which requires almost twice the amount of 
uranium than the mini subcritical reactor. 
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1. Introduction 

Nuclear energy is a reliable option for worldwide 
energy supply, also is a good option to mitigate the 
greenhouse gas emissions; however, it is perceived as 
risky energy source [1]. 

Some developed countries like Japan, Germany, 
Switzerland, etc. are limiting their expansion on nuclear 
power production. On the other hand, countries like India and 
China are increasing the nuclear power plants construction 
due to the energy demand combined with the shortage of 
fossil fuels and the limits on greenhouse gas emissions [2]. 
In the case of countries in Latin America, nuclear energy 
production is very small, however expectations of increase 
its participation are promising. 

Most countries meet their energy demands based on 
fossil fuel consumption, contributing to global warming 
through the emission of gases such as sulfur dioxide, 
nitrogen oxides and carbon dioxide. To avoid this has 


been used the sources of alternative energy such as wind, 
hydraulic, solar and nuclear energy etc. But none of these, 
in isolation, can completely satisfy the energy demand. 
Worldwide there are 447 nuclear power plants in operation 
and 287 new will be on line by 2035 [3, 4]. 

There are many ways to classify nuclear reactors, in 
terms of their application nuclear reactors are classified as 

[5], 

• Power reactors that are used to produce electric energy, 
through the nuclear fission chain reaction induced 
by neutrons in a multiplied medium. Commercially 
reactors are designed to produce electric power of 600 
and 1500 MWe with thermodynamic efficiencies close 
to 33% [5]. 

• Research reactors are used in several applications and 
to carry on scientific research. This type of reactors 
becomes a source of neutrons that are used in the 
analysis of materials by neutron activation analysis or 
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to produce radionuclides being used in medicine 
and in industry. The power of these reactors varies 
from a few Watts to almost 100 MW [6-8]. 

• Teaching reactors, which are normally subcritical 
meaning that the nuclear fission chain reaction 
occurs only with an external neutron sources; as 
soon this source is withdrawn or it is switch off the 
reactor is shutdown. This type of reactor are easy 
and safe to operate [7-10]. 

Subcritical reactors are used to perform experiments that 
confirm calculations of power reactors, and are also being 
used to transmute some radionuclides in nuclear waste 
[10, 11]. The most important parameters of a reactor are: 
the multiplication factor (k ), the reactivity, the spatial 
distribution of neutrons, the neutron spectra inside and 
outside the reactor core and the ambient dose equivalent 
outside the reactor [12,13]. 

The k eff is the ratio between the neutron density pre¬ 
sent in the i th + 1 generation and the neutron density present 
in the i th generation. If k gff > 1, the reactor is supercritical, 
if, k gff = 1 the reactor is critical, and if k eff < 1 the reactor is 
subcritical [14]. Subcritical reactors are designed to remain 
subcritical regardless of any situation like flooding or if the 
external neutron source is left inside [15]. 

In Mexico there are two subcritical assemblies, one 
in the Instituto Politecnico Nacional (IPN) in Mexico 
City and the second one is in the Universidad Autonoma 
de Zacatecas (UAZ). Both reactors are Nuclear Chicago 


model 9000 which is a heterogeneous reactor using light 
water as moderator, and natural uranium as fuel. Moderator 
and fuel are contained into a 121.92 cm-diameter and 
152.4 cm-height stainless steel cylinder. Nuclear fuel 
is distributed into a hexagonal array in whose center is 
inserted an isotopic neutron source to start-on the reactor. 
The neutron source is a 1.01E(7) s 1 239 PuBe [16]. With 
all fuel elements (2.43 tons of uranium) the k eff is 0.8268, 
having an amplification factor p = 5.77, the reactor power 
is 0.617 mW, the maximum neutron flux is cp M . x = 2.93E(4) 
cm -2 -s '. The ambient dose equivalent outside the stainless 
steel container is 14 microSv/h [7]. 

The design of nuclear reactor can be carried out using 
deterministic or Monte Carlo methods. The deterministic 
methods are used for calculations of critical reactors, these 
methods are approximations of the Boltzmann equation, 
like the multi-group diffusion theory [17-19]. 

The objective of this work was to use Monte Carlo 
methods in order to design a subcritical mini nuclear reactor 
with natural uranium fuel, moderator of light water and a 
neutron isotope source of 239 Pu-Be, and to compare this 
design with the Nuclear Chicago, model 9000, subcritical 
reactor. 

2. Materials and Methods 

This work was carrying out with the Monte Carlo code 
MCNP5 and the cross section library ENDF/VI [20]. The 
reactor was designed modeling the nuclear fuel concentric 
rings, as is shown in Figure 1. 


Fuel Moderator 

[Wi 

Neutron source 


Figure 1: Side view of nuclear reactor model with five rings of uranium. 


Figure 1 shows the reactor model with five fuel rings, 
in the center is the 239 PuBe source, and the moderator is 
light water. Five cases were designed having 8.5, 14.5, 
20.5 26.5, and 32.5 cm-inner radius, each fuel ring is 3 
cm-thick, and 36 cm-high. 

The Monte Carlo method was used to calculate 
the effective multiplication factor, k eff5 that was used 
to calculate the neutron amplification factor, p, using 
equation 1. 

1 

// =- 

(i) 


Also the spatial distribution of the neutrons and the 
ambient dose equivalent was estimated, as well as the 
neutron spectra along the radius for reactor with one, 
two, three, four, and five fuel rings. With the k eff and 
the p values the reactor power was calculated using 
equation 2. 


V p f 

In equation (2), Q is the neutron source strength 
1.01E(7) n/s, v = 2.44 n/fission, andp f = 3.20077E(10) 
fissions/Watt. 
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Monte Carlo calculations were also used to 
estimate the neutron spectra to different distances 
in the reactor for the 239 PuBe source, and the reac¬ 
tor with 1, 2, 3, 4, and 5 rings of fuel (RF). 

The ambient dose equivalent rate was 
calculated to 51 cm using the neutron fluence- 
to-ambient dose equivalent coefficients from the 
ICRP 74 [21]. 


3. Results and Discussion 

In table 1 are shown the values of fuel mass, k cff , amplifica¬ 
tion factor (p), reactor power, and the ambient dose equiva¬ 
lent rate H*(10) in function of the amount of reactor fuel 
rings and when in the reactor center is inserted a 239 PuBe 
source emitting 1.01E(7) n/s. The H*(10) values are to 51 
cm from the center (this is 1 cm outside the reactor surface). 


Table 1: k eff , fuel mass, p, reactor power and H*(10) for reactor with different fuel rings. 


Rings of 
fuel 

Fuel mass 
[kg] 

keff 

F 

Power 

[mW] 

H*(10)(51 cm) 
[jxSv/h] 

1 

90.7 

0.5298±3.494E(-4) 

2.13 

0.27 

13.4 ±4.2% 

2 

259.2 

0.6530±6.610E(-4) 

2.79 

0.43 

14.4 ± 5 . 7 % 

3 

505.5 

0.7186±5.617E(-4) 

3.70 

0.63 

18.3 ±5.3% 

4 

829.5 

0.7585±5.648E(-4) 

4.42 

0.79 

29.4 ±3.9% 

5 

1229.2 

0.7819±5.763E(-4) 

5.01 

0.93 

53.0 ±3.25% 


The k ff , p and the reactor power increases as 
the amount of fuel is added. The increase of k „ 

eff 

with fuel mass has the same shape reported for a 
Nuclear Chicago subcritical reactor. In the case 
with 5 rings of fuel the k eff is 0.7819 the total 
mass of uranium is 1229 kg, this reactor is a 50 
cm-radius and 50 cm-height cylinder. 

The k of a Nuclear Chicago subcritical 
reactor is 0.8268 which is 5.7% larger than the 
k „ of our mini reactor. This difference is small 

eff 

considering that the Nuclear Chicago subcritical 
reactor requires 2430 kg of fuel that is inside a 
61 cm-radius and 152 cm-height cylinder. The 
amount of fuel of Nuclear Chicago subcritical 
reactor is almost twice (1.977) the amount of 
fuel of mini subcritical reactor. 

The volume of Nuclear Chicago subcritical 
reactor is 4.525 times larger than the mini 
subcritical reactor. 

The amplification factor of Nuclear Chicago 
subcritical reactor is 5.77 being 15.2% larger 
than p in the mini subcritical reactor. On the 
other hand the reactor power of mini subcritical 
reactor is 0.93 mW which is 50.7 % larger that 
the power of Nuclear Chicago subcritical reactor. 

Outside the reactor (1 cm form the external 
surface of reactor) the ambient dose equivalent 
rate in the mini subcritical reactor is 53 mi- 
croSv/h which is 3.8 times larger than the H*(10) 
in the Nuclear Chicago subcritical reactor; this 
difference is probably due to the reactor size. 

In figures 2, 3, 4, 5, and 6 are shown the 


neutron spectra to 3.8, 10, 20, 30, and 40 cm 
from the neutron source center along the radius 
respectively. Neutron spectra are for cases where 
there is only the neutron source (no fuel), and 
when the rings of fuel (RF) are added. 

In figure 2 are the neutron spectra to 3.8 cm- 
radius and in the middle point along z-axis. The 
bare 239 PuBe isotopic neutron source produce 
neutrons from 0.5 up to 11.5 MeV, to 3.8 cm in¬ 
side a light water moderator epithermal and ther¬ 
mal neutrons are produced due to the moderation 
and thermalization process. To 3.8 cm the neu¬ 
tron spectra have the same shape but the amount 
of neutrons increases as the fuel rings are added. 
The spectra have a small peak of thermal neu¬ 
trons (E < 4 eV) and a larger peak of fast neu¬ 
trons (0.1 < E < 20 MeV) linked by epithermal 
neutrons (4E(-7) < E < 0.1 MeV). 

In figure 3 are the neutron spectra to 10 cm 
from the reactor center, which is the moderator 
located between the first and the second fuel 
ring. 

To 10 cm the impact of neutron multiplication 
due to nuclear fissions is larger than at 3.8 cm. 
Bare source has the largest amount of thermal 
neutrons in comparison to spectra produced when 
ring fuels are added; this is due to the capture 
of thermal neutrons in the fuel. As the RF are 
included the fast neutrons are shifted to the left 
where the larger peak is for 0.5 to 1 MeV, while 
the bare source has its largest peak in 2 to 3.5 
MeV neutrons. 
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Figure 2: Neutron spectra, per history, to 3.8 cm from the center 



Figure 3: Neutron spectra, per history, to 10 cm from the center 


In figure 4 the neutron spectra are to 20 cm-radius this 
is inside the third fuel ring. To 20 cm neutron spectrum 
due to bare 239 PuBe is similar to the neutron spectrum 
when the reactor has 1 RF for neutrons whose energy 


is larger to 0.4 eV, but thermal neutron group is larger 
when the reactor has 1 RF. With 2 RF the amount of 
thermal neutrons is larger in comparison to reactor with 
1, 3, 4, and 5 RF. 



Figure 4: Neutron spectra, per history, to 20 cm from the center 
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Neutron spectra to 30 cm, which is located inside the but when the fifth RF is included, in this case the larg- 
fifth fuel ring, are shown in figure 5. In this point the est contribution is due fast neutrons probably due to 
amount of thermal neutrons is the largest for all cases the amount of neutrons produced during fission. 



Figure 5: Neutron spectra, per history, to 30 cm from the center 


The fifth RF has a radius of 32.5 cm, and the reactor 
radius is 50 cm, the space from 32.5 cm to 50 cm is 
filled with water. Neutron spectra to 40 cm, which is 
in the moderator, are shown in figure 6. In this point 
the main contribution are thermal neutrons where the 
amount is increased as the fuel rings are added. The 
reduction of fast neutrons and the increase of thermal 
neutrons as the distance is increased are similar to 
measurements carried out in a Nuclear Chicago model 
9000 subcritical reactor [22], and a subcritical reactor 


with TRIGA reactor fuel [23]. 

In the middle point of reactor height total neutron 
fluence was calculated, in figure 7 is shown the neutron 
distribution when the reactor has the neutron source 
and the first fuel ring. 

In figures 8,9,10, and 12 is the neutron distribution 
for 2, 3, 4, and 5 RF reactor. In all cases the amount 
of neutron is small outside the reactor (51 cm). The 
neutron distribution is strongly affected by the amount 
of fuel rings in the reactor. 



10* 10' 10* 10* lO^ 1 10* 10* 10 1 W 10 1 10 ? 


Energy [ MeV J 

Figure 6: Neutron spectra, per history, to 40 cm from the center 
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Figure 7: Neutron fluence distribution for reactor with 1 RF 



Neutron fluence 
per history 
[ n/cm 2 ] 


Figure 8: Neutron fluence distribution for reactor with 2 RF 



Figure 9: Neutron fluence distribution for reactor with 3 RF 
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Figure 10: Neutron fluence distribution for reactor with 4 RF 



Figure 11: Neutron fluence distribution for reactor with 5 RF. 


In figures 7, 8, 9, 10 and 11 the distribution of neutron 
fluence per history are shown for mini subcritical 
reactor with 1, 2, 3, 4, and 5 fuel rings, the shape of 
neutron fluence distribution is alike to the distribution 
reported for the Jordan subcritical reactor [24]. As 
the rings of fuel are added the neutron distribution is 
enhanced due to neutrons produced during fission. In 
any point in the reactor the neutron fluence is the results 
of neutrons produced by fission that are transported 
less those neutrons that leak out from the reactor and 
those absorbed by the fuel and the moderator. 

4. Conclusions 

Monte Carlo methods were used to design a mini 
subcritical nuclear reactor. This reactor has a 239 PuBe 
isotopic neutron source and has light water as 
moderator. The reactor is heterogeneous where the fuel 
is cylindrical rings. In the design five reactors were 
calculated using 1, 2, 3, 4, and 5 fuel rings. The k eff , 


jii, the reactor power and the ambient dose equivalent 
rate were calculated, as well as the neutron spectra and 
the neutron distribution. The features of reactor with 5 
rings of fuel were compared with a Nuclear Chicago 
subcritical nuclear reactor. 

All mini subcritical nuclear reactors have a k „ less 

eff 

than 1, ranging from 0.5298 to 0.7819, the increase of 
k is due the increase of fuel mass [25]. 

Neutron spectra have thermal, epithermal and 
fast neutrons. The amount of neutrons for each energy 
group depends upon the site where the spectrum is 
calculated. 

Near the neutron source the amount of fast 
neutrons is larger than thermal neutrons, as the 
distance is increased the amount of thermal neutrons 
is increased. 

Neutron fluence distribution tends to enhance as 
the amount of rings of fuel are added. 

Considering the volume and fuel mass of mini 
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subcritical nuclear reactor with 5 rings of fuel in 

comparison to volume and fuel mass of a Nuclear 

Chicago subcritical reactor the difference of k eff is 

small. 
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